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We formulate a model for the steady state response of a nonlinear quantum oscillator structure, 
such as those used in a variety of superconducting qubit experiments, when excited by a steady, 
but not necessarily small, ac tone. We show that this model can be derived directly from a circuit 
description of some recent qubit experiments in which the state of the qubit is read out directly, 
without a SQUID magnetometer. The excitation profile has a rich structure depending on the 
detuning of the tone from the small-signal resonant frequency, on the degree of damping, and on the 
excitation amplitude. We explore two regions in detail: First, at high damping there is a trough in 
the excitation response as a function of detuning, near where the classical Duffing bifurcation occurs. 
This trough has been understood as a classical interference between two metastable responses with 
opposite phase. We use Wigner function studies to show that while this picture is roughly correct, 
there are also more quantum mechanical aspects to this feature. Second, at low damping we study 
the emergence of sharp, discrete spectral features from a continuum response. We show that these 
the structures, associated with discrete transitions between different excited-state eigenstates of the 
oscillator, provide an interesting example of a quantum Fano resonance. The trough in the Fano 
response evolves continuously from the "classical" trough at high damping. 



I. INTRODUCTION 

Spectroscopic studies of superconducting qubit systems have increased steadily in sophistication in recent years. For 
many of the first Josephson-junction qubits, their only link with the outside world was magnetometry: diamagnetic 
currents, particularly in flux qubits, were detectable by SQUID magnetometers. While this provided the basis for 
some of the first experiments in which solid state qubit states were detected, these magnetometers were very intrusive 
devices in these systems, and were clearly an impediment to scalable systems. Spectroscopy with these SQUID 
detectors was very limited. 

With the application in recent years of many sophisticated techniques from microwave engineering, various alterna- 
tives to the SQUID approach have emerged. The detection, protection, and manipulation of superconducting qubits 
have undergone impressive improvements. In particular, the spectroscopic probing of qubits and qubit /resonator 
systems is now capable of very useful high-fidelity qubit detection)]] [SJ, as well as detailed characterization of the 
parameters of these systems opening a beautiful view [3] into the rich quantum response that these systems exhibit. 

This paper will provide a theoretical study of a model that is central to these spectroscopic studies. The model is 
that of a single oscillator (bosonic mode) with a purely quartic anharmonicity. The oscillator is driven by a single ac 
tone, which is not necessarily a small perturbation on the system, and the system is damped with a purely frictional 
term, representative of a zero-temperature quantum environment. The observable is taken to be essentially the 
displacement of this resonator. 

We will emphasize here the insights that can be gained from analytic or semi-analytic studies of this model. In 
previous studies of this model[l], it has been shown that the steady state response of the system can be calculated 
in closed form. It is very rewarding to look at the many aspects of this amazingly rich solution. At high damping it 
exhibits much of the phenomenology of the classical Duffing oscillator. In a broad region of detuning, on the "red" side 
of the harmonic response frequency, and for moderately high amplitude excitation, the response is "high" , or quasi- 
resonant. But as damping is decreased, new, sharp resonant structures emerge from the broad Duffing feature. These 
lines are clearly associated with individual spectroscopic transitions in the anharmonic ladder of energy eigenstates. 
But the line shapes of these new lines are very non-Lorentzian, and are indicative of a complex interaction among 
different response channels. 
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In the present work we do further calculations to gain more understanding of several of these spectroscopic features. 
Up to a point, the behavior in the high damping regime exactly matches the classical Duffing theory. There is a linear 
response regime, in which the response is a simple Lorentzian centered on the harmonic frequency of the oscillator; 
at higher excitation amplitude the resonant line moves to lower frequency and sharpens. But there are differences: 
the quantum model intrinsically does not have multiple solutions for the steady state response, and so does not have 
any hysteretic behavior. Near the classical bifurcation point the quantum response shows a fairly abrupt transition, 
as a function of excitation detuning, between a "low" and "high" response regime. 

While some details of this interesting quantum behavior have been discussed in previous work, we shed new light 
on this behavior using simple Wigner function computations. In the middle of this transition from "low" to "high" 
response, there is a dip in the response amplitude, which has been interpreted [3J H] as a kind of intermittency in 
the steady state involving two responses with different, and therefore canceling, phases. Our Wigner functions show 
that there is clearly some truth to this point of view, but we find that, even for damping large enough to wash out 
all other quantum spectral features, the dip has more features than can be explained by any classical intermittency 
point of view. We find that this story is considerably enriched by using some of the modern tools of quantum 
information theory, such as calculations of the von-Neumann entropy, and some basic facts about the eigenanalysis 
of the completely positive map describing time evolution in our model. 

We also explore in detail how structure emerges as the damping parameter is decreased. While the spectral response 
grows considerably in complexity, we find that the simple "dip" described above in the high damping regime persists 
for all values of damping, but comes in a different guise in this quantum limit. In one part of the spectrum, we identify 
this dip as a feature of a Fano resonance. The Fano phenomenon is a quintessentially quantum interference, yet we 
find that in this case it is one extreme of a family of interference phenomena that also include a nearly classical form 
of interference at the other extreme. 

These studies help us to answer a very fundamental question about anharmonic resonant structures, which is: when 
is it reasonable to call such a structure a qubit? It is evident that the answer is a dynamical one: the very same 
resonator can be either a very well defined two-level quantum system (a qubit), or an oscillator requiring a many-level 
quantum description, depending on the amplitude and frequency of oscillation. In this paper we will fairly precisely 
demarcate the rather complex boundary between these two regimes. 



II. THE ELECTRICAL DEVICE 
A. Circuit basics 

The calculations here are motivated by current experiments on superconducting circuits such as the one shown in 
Fig. [lj The ac voltage signal (in the microwave frequency range, in cases of current interest) excites a transmission 
line; this transmission line is coupled capacitively to a small, discrete superconducting circuit that we will, for short, 
call "the qubit." But an important message of this paper is that structures of interest do not always behave as qubits, 
that is, like two-level quantum mechanical systems. From the point of view of electrical science, the combination shown 
of a capacitor, inductor, and Josephson junction is viewed as an anharmonic resonator. As we will show, depending 
on the excitation conditions, this classical view is sometimes more appropriate than the quantum mechanical one. In 
fact, there is a rich set of phenomena that are explained by neither the classical or the simple (two-level) quantum 
point of view. 

Many features of a real experiment are omitted from Fig. [T] Naturally the qubit is held at low temperature, while 
the microwave voltage source is (presently) always at room temperature; the diagram does not show the attenuator 
stages that are necessary to keep room temperature noise out of the qubit. But these stages are engineered so that 
the functioning is that of a cold, matched (i.e., Zg — Zq) source. On the output end, conceptually measurements of 
the qubit are recorded as the voltage across the terminating resistor Z^. In reality, this circuit is a complex, active 
chain of amplifiers and other components, with a signal recorded at room temperature. But again, the engineering 
of this chain has the object of reducing the effective functioning to that of the simple circuit in Fig. [T] in which the 
load is matched (Zl = Zq) and cold. 

Much experiment and analysis (c/. [S]) has been devoted to a related circuit, in which the transmission line is 
only weakly (capacitively) coupled to the source and load circuits. This is the "circuit QED" setting, in which 
the transmission line forms a resonator which functions as a quantum-coherent resonator. The combination of this 
quantum resonator and the qubit exhibits rich, interesting, and potentially useful physics. The system we analyze here 
does not have all the same potentialities as the cQED system, but its dynamics still has a great deal of complexity in 
it, and the simpler circuit has its own potential application in the realm of simple, reliable qubit measuring circuits. 




FIG. 2: Circuit for experiment. Zq = 5011. 
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As far as the dynamics of the qubit is concerned, the circuit of Fig. [2] accurately describes the schematic of Fig. 
[I] in the matched condition Zq = Zg = . The voltage V% at the node indicated is the observable of the actual 
experimental setup. We have replaced the inductor and Josephson junction by a nonlinear inductor, which is sufficient 
given that we will not wish to describe the global response of the device (e.g., to large changes of an external bias 
flux). We have also included resistor R, to describe internal loss of the qubit. This is in some cases competitive with 
the losses due to the "external" resistors in the circuit of Fig. [2] 

We can get another circuit more directly amenable to a quantum treatment by using some circuit transformations 
due to Zmuidzinas et al.[6 , and making a "narrow band" approximation, which amounts to saying that the only 
response of interest in the Fourier transform of the voltage will be at frequencies near the excitation frequency; this 
will be compatible with the rotating wave approximation that we will introduce later. This is shown in Fig. [3j 




V 



FIG. 3: Equivalent circuit to Fig. [2] Note that R e 
then identical in the two circuits. 



2Z /(uj CcZo) 2 with cjq = 1/y/LC and I eq = V s /Z . The voltage V is 



Note that the measurable quantity V2 has disappeared from this circuit. But is easily recovered from the voltage 
V via the linear circuit relation[3] 

v 2 = \v s + (iu C c Z )V (1) 

This is valid for loqC c Zq « 1. Eq. ([I]) is understood to be in the frequency domain, with the "?" indicating that the 
contribution of V to the signal V2 is out of phase with that of the drive V s . 

We will take V s to be a classical ac variable with a simple sinusoidal time dependence: 

V r a (t) = V r »coB(wp*) (2) 

We will assume that the nonlinear inductance is characterized by the small-signal energy storage formula 

1$ 2 1$ 4 

^2T + 4Ll ® 

$ is the (formal) magnetic flux through this inductor, i.e., $(t) = J* Vi,(t')dt' . Other powers of $ will generally 
appear in this energy, but our truncation will often be very accurate. This energy expression implies the non-linear 
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inductive two-terminal response 



/l - a* - l + u ■ (4) 



Using any standard circuit mechanics treatment (e.g., [7]), the Hamiltonian of the circuit of Fig. [3] is then 

2 C" V 2 £ 4L 4y / Z W 
Here C" = C + C c . The classical equation of motion for $ will also have a "friction" term going like 



$ 1 1 1 2Z Q 

_i_ ^ 

i?'' R' R eq R 7 eq (ojqCcZq) 



Note that the voltage V of interest here is related to our canonical momentum via the two-terminal relation for the 
linear capacitive element: 

V-±. (7) 
B. Quantum treatment 

We go over now to a quantum treatment of the circuit. This will give Q, and therefore V and V2, the status of a 
quantum operators. We will assume that the experiment consists of a sequence of weak measurements of V2, so that 
we will extract the expectation value (V2), and that the system is not otherwise disturbed by the measurement. The 
standard quantum commutation relations are 

[*,Q]=»fi (8) 
Introducing bosonic creation and annihilation operators, we have 

and similarly for 4>. Here the subscript on the operator a stands for "Schrodinger picture," as we will soon do most 
of our calculations in an interaction picture. 

Using this operator expression for Q in Eqs. ([9]), we can get a quantum expression for the observable V%. Going to 
the time domain: 

(V2) = \v s cosujpt + CJ0 ^^ (Re(a) cosw p t + Im(a) sin (10) 
The Hamiltonian operator describing the coherent part of the system evolution becomes 



H s = hu al ch a Sch + — (a Sc h + ag ch ) + 2Hecos(ui p t)(a Sc h + ag ch ) (11) 
We now go to an interaction picture (a = e _,Wp ag^, etc.), and perform a rotating-wave approximation, obtaining 



ff = /iAo t a + lix(o t ) 2 a 2 +l!e(o + ffl t ) (12) 
Here the detuning A is A = luq — uj p , and we introduce scaled parameters, all with units sec -1 : 

"°-7m' x -tc^> e -vm[v) y (13) 
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The lossy parts of our circuit are accounted for using the usual bosonic bath description; making the standard Born- 
Markov approximations OH], and assuming the bath to be at zero temperature, the system dynamics is described 
by the Lindblad master equation: 

g = -i [H, p] + 1 (2apJ - a) ap - pa) a) (14) 

with 

^ = ^Cr < 1B > 

In this rescaled form, our model is identical to one employed recently in nanomechanics by Babourina-Brooks, Doherty, 
and MilburnjS]. In fact, its essential features were already calculated thirty years ago by Drummond and Walls [1]. In 
this paper, we will explore this rich solution and provide physical insights into some of its many interesting regimes. 

C. spectroscopy 

We will model the "spectroscopy" experiment as done in [1] and in many other recent studies of superconducting 
qubit systems. Spectroscopy is measured as the magnitude of the ac output voltage signal at the same frequency as 
the input. For the model above, this will be proportional to the expectation value of a in steady state: 

(a) = Tr( a p ) (16) 

Where po is the steady state response: 



dpo 
dt 



(17) 



We can write the Lindblad master equation in the form 

where S is a linear superoperator, generating a TCP (trace preserving completely positive) map e St . The resulting 
equation for the steady state, 

S\po] = (19) 

states that po is obtained as the zero eigenvector of the linear non-Hermitian (super)operator S. It is this point of 
view that is employed in the numerical studies that we present below. 

But numerics are only necessary for developing physical explanations, since the spectroscopy calculation itself, 
remarkably, has a closed form solution. Quoting Eq. (5.17) of [3]: 



0F2 



A+x-n/2 A-N7/2 2£ 
X ' X ' X 2 



A - l"f/2 j? ( A-ij/2 A+i-,/2 2£ 



X ' X ' X 



(20) 



Here 0-F2 is a standard hypergeometric function. 



III. RESULTS 



For three different damping parameters 7, Fig. [4] shows \{a}\, which is proportional to the magnitude of the 
transmitted voltage in a spectroscopy experiment. The plots are shown as a function of A, the detuning away from 
the harmonic resonance frequency, and e, the drive intensity. These quantities are normalized to the anharmonicity 



parameter x (i.e., X — 1 m Eq. (12)). The three panels show three different settings of the damping parameter 



(again, relative to %). For the largest 7 shown (7 = 2 in part (a)) the response is in some sense classical; no structure 
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FIG. 4: Response amplitude |(a)| as a function of detuning A and excitation amplitude e, all for anharmonicity parameter 
X = 1. (a) Moderately high damping (7 = 2). Red line: scan across the response at e = 3.2. Points A, B, C, and D along 
this scan path are examined further in Fig. [5] The trough at point B is easily visible in the line scan. Dark lines show 
the boundary of the metastable points [10 for the classical Duffing oscillation with these same parameters, (b) Response for 
moderate damping 7 = 0.3. Trough has become more prominent but is now scalloped, indicating the development of discrete 
quantum resonances, (c) Response for low damping 7 = 0.01. Quantum resonant structure becomes very prominent, (d) A 
closer view of the first resonant structure in the low damping case 7 = 0.01. 

involving transitions between individual quantum mechanical levels is present at any excitation level. As we will 
discuss, though, there are many aspects of this result that are not classical in any naive sense. 

The principal feature of Fig. |4|a) at moderate excitation levels (e > 1) is a fairly abrupt upward step in the 
response, with a brief dip in between. This dip or trough has been noted previously [2H1]j an d is prominently seen in 
recent experiments; we will discuss its origin shortly. As 7 is decreased (Fig. (4^b)), this trough develops ripples that 
are synchronized with the quantum transitions in the anharmonic ladder of states. At very small 7 (Fig. |4|c)) these 
ripples resolve themselves into spectroscopic resonant peaks, which merge into the "high" continuum at large e. As 
we will discuss later, it is only in this third setting that it is meaningful to say that the anharmonic resonator behaves 
as a qubit, for the most part only at small excitation levels e < \ and near zero detuning |A| < \- 



IV. RESULTS: HIGH DAMPING REGIME 



The true classical behavior would be quite different from what is seen in Fig. Qa), as was noted already in (3). 
Inside the two lines shown in Fig. |4|b) , these are two classical steady state solutions - the dynamics bifurcates at the 
meeting of these two lines, and is hysteretic above the bifurcation. A quantum treatment, even one with large loss, 
cannot show this bifurcation, although the quantum solution is in fact very close to the classical solution in the "low" 
region beneath the lower bifurcation line, and in the "high" region above the upper bifurcation[4 . Qualitatively, this 
difference can be ascribed to the possibility of transitions between the two classical steady states; mathematically, 
it arises because the superoperator S[p] has a unique zero eigenvalue. Classical hysteretic behavior can be and is 
recovered in a quantum treatment of the transient behavior, as studied recently in 



Collapse of hysteresis due to quantum effects is actually well known [T2]. Special to the anharmonic resonator 
system, however, is the robust trough easily seen in Fig. |4|a). As originally noted in [4], the trough can be reasonably 
explained by the observation that the classical "high" and "low" states have almost opposite phases, since the first is 
essentially a resonant response, while the second is non-resonant. 

In the simplest view of the trough, the state intermittently switches back and forth between the high and low 
Gaussian state. While this gives a rough description of the situation at point B in Fig. Qa), the state there is in 
fact not purely a mixture of a low-amplitude with a high-amplitude Gaussian state. We explore this transition by 
examining the Wigner function W of the steady state p as we pass through the "trough" by varying the detuning 
A. We use the standard formula 14J 

W(a) = Tr (D(-a)p D(a)Il) 

OC 

D{a) = exp (oa* - a*a) , II = ^(-l) n |n)(n| (21) 

with complex displacement a. The simplest is the Wigner function at the "low response" point A (Fig. ga)). To the 
eye, this state is nearly a coherent state slightly displaced (i.e., with substantially less than one photon of excitation). 
In fact, the state has some excess noise (i.e., has a somewhat broader Wigner function than a coherent state), which 
can be measured by the von Neumann entropy Svn(p) = — Trplogp. At point A, 5 i v n(po) — 0-03 bits. 




FIG. 5: Wigner functions for four points shown in Fig. [4] 

The Wigner function at the "high response" point D departs further from the simple classical view. This Wigner 
function (Fig. [5](d)) of course has a larger \(a)\ than in point A, but it is much farther from a coherent state. Despite 
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its appearance, it is not number squeezed (the radial width is not smaller than a coherent state). It is perhaps better 
described as a coherent state with significant excess phase noise. The entropy of this state is 5 v n = 0.85 bits. 

The Wigner function in the transition region deserves considerable discussion. We will focus in detail on point C. 
The classical view that the state is a mixture of two classical responses is not entirely incorrect. But it fails to capture 
many of the aspects of the actual response. The system clearly spends appreciable time at phases intermediate between 
the high and low response, with n = 2.56, so it should not be viewed as quickly switching between the two. The state 
has high entropy, as a classical mixed state would, but its precise value (SVn = 1-74) gives us more information about 
the nature of the state, as we explain now. 
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FIG. 6: Wigner functions of the extremal metastable responses at point C (A = —5.2) in the response diagram, showing that 
the picture that the Wigner function of the response (Fig. [5jc)) involves two possible metastable states is approximately, but 
not precisely, correct. Color scale is the same as in Fig. [5] 




FIG. 7: von Neumann entropies for the steady state p and the extremal me tastable states p+ and p_, as a function of A for 
e = 3.2. Also shown are the coefficients f3±, suitably normalized, in Eq. (23 1. 
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The mixing picture is illuminated by a computation which pulls out a version of the two "classical response" states 
from the full state of Fig. [5^c); this is shown in Fig. [6] A vestige of the bistability of the classical response is clearly 
reflected in the eigenstructure of the superoperator S at the operating point C. The second smallest eigenvalue if 
S is at Ai = —0.215, while the third eigenvalue is at A2 = —2.204. It is the presence of two very nearby low lying 
eigenvalues Ao,i, with a large separation from the rest of the spectrum, is a reflection of the classical bistability. It 
is straightforward to calculate this second eigenvector |1) a.k.a. 8p\: S\l) — Ai|l). We use the notation Sp as a 
reminder that while this eigen-matrix is Hermitian (because it is associated with a real eigenvalue, rather than with 
a complex-conjugate pair[]l)]), it is traceless. But the important object we can construct from Spi is 

Po+PSft. (22) 

For real (3. There will be a range of /3s, including (3 = 0, for which the operator will be positive semidcfmite. We find 
that the two extremal operators, the one with the smallest possible j3 = (3- < and the largest possible (3 = /3+£0, 
are the ones corresponding to the two classical solutions. Fig. [6] shows the Wigner functions of these two density 
operators, 

P±=Po + f3±S P i (23) 
We can see that this does a reasonably good job of separating the response into a high-amplitude and a low-amplitude 



part. This is also confirmed by the calculation of the von Neumann entropy of the range of states Eq. (22) (Fig. 
W\. The two extremal states p± are those of the lowest entropy. Classically, the states in the interior of this range 
(including po) should have an excess entropy above the weighted average of the extremal states (see Fig. [7]) equal to 
the ordinary binary entropy of mixing H(x) = —a; log a; — (1 — x) log(l — x). In fact Fig. 7] differs from this classical 
situation in two respects: 1) The maximal entropy of mixing is not 1 bit, but about 0.74 bits. 2) The excess entropy 
is nearly proportional to H(x), but not exactly so, reflecting a quantum effect, the noncommutivity of the extremal 
states p + ] 7^ 0. 



V. RESULTS: LOW DAMPING REGIME 



We now turn our attention to the lightly damped regime (by which we will mean 7 < We expect that it is in 
this regime that our anharmonic oscillator "is" or "behaves like" a qubit. We will explore the light that our response 
calculation can shed on the question of when such an oscillator is or is not a qubit. 

The outstanding qualitative feature of the response as 7 is lowered (as in the series in Fig. [4]) is that stalactites 
form, descend, and sharpen from the bottom of the high-response continuum. The continuum itself, at least in the 
vicinity of the classical upper bifurcation line of Fig. |4]Ja), is not much affected by the lowering of 7. The trough 
remains an important feature, becoming first slightly, and then strongly, scalloped, finally breaking up into segments 
between one stalactite and the next. 

Given that the quantum energy levels to n for the quartic oscillator obey to good approximation the simple rule 

uj n+1 - u n = cj + nx, (24) 

it is natural to associate the resonant structure at zero detuning with the "qubit" — > 1 transition, the structure at 
detuning A = — x with the 1 — > 2 transition, that at A = —2% with the 2 — > 3 transition, and so forth. We might 
be tempted to guess that excitation detuning A = — nx, the system is acting like a qubit involving the levels \n) and 
\n + 1). As we will seen this is most definitely not the case, except for n = 0. A detailed inspection of the n = 
"stalactite" shows that its structure is quite distinct from that of all the others. In fact, the excitation structures at 
higher n are indicative of interesting multi-level behavior, as we will see. 



A. A perturbation theory 



We will do a perturbation analysis in e to elucidate the nature of the resonant structures that appear as the 
excitation level is increased. While much can be learned simply by Taylor expansion of the exact expression for |(a)|, 
much more is revealed by developing perturbative expressions for the eigenoperator of the superoperator S. 
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To develop this, we first record an explicit matrix form for this superoperator: 

Ski,ijpij = {-i[A(k - I) + X [k(k - 1) - 1(1 - 

-ie(Vk + 16k t i-iSij + Vk5k,i+i5ij - VlS kl S j+ ij - Vl + I5ki§j-i,i) 

+ |[2 v /(fc + l)G + l)4+i^/+i J - (k + 04;<%] }py (25) 

S is sparse, but of course it is not, as a matrix, Hermitian. We analyze the steady state response by a perturbation 
approach, we can separate S thus: 

S = S + eV (26) 

Sq can be analytically diagonalized - it has a "bidiagonal" form as a matrix. It is convenient to index its eigensolutions 
by two integers n and q. Since So is non-Hermitian, it is necessary to distinguish its right from its left eigenvectors, 

So\nqn) = Kiq\nq R ) (27) 



{nq L \S = (nq L \\ nq (28) 

Recall that the left eigenvectors are not the adjoint of the right eigenvectors, but the matrix formed by the left 
eigenvectors is the inverse of the matrix of right eigenvectors. The physical steady state response to zeroth order in e 
is the eigenvector |00_r). 

So has both real and complex eigenvalues, for both of which ReA ng < for e ^ (except for Aoo which is always 
equal to zero). It is straightforward to confirm the following: All the complex eigenvalues of So may be written 

Kg = - (q + |) 7 - niA - n(n - 1 + 2q)i x (29) 

with integers q > 0, n > 0. The right eigenvectors \nq) (a.k.a. p nq , see Ref. [T3]) are operators on the Hilbert space 
with matrix elements 

(-l-2n«Y 

{s\p„ q \t) = 6 s , t+ ± — . (30) 

(q — tjWs'.t'. 

Each has a complex conjugate partner \ nq = A* g , p nq — p^. The real eigenvalues and right eigenvectors are given 
by the same formulas with n — 0. We do not record the left eigenvectors (nq\ here, but they are as straightforward 
to write down as the right eigenvectors. 

The right eigenvector |00^J n ) of the full superoperator can be conveniently developed using Brillouin-Wigner per- 
turbation theory [To]. Going over to a schematic single-integer indexing of the eigenstates: 

\0 R ) = |0 R ) + e^ +e ^ (Eo-XMEo-Xj) + ~ (31) 

fe>o fe>oi>o v u M ■>' 

The Brillouin-Wigner formula is frequently not used because the exact eigenvalue Eo is unknown. But in this 
application, of course, E = exactly; so this series is actually very useful, and its first few terms are readily 
calculated. 



B. Characteristics of the resonant response 



Using this perturbation theory, we can write out in matrix representation (in the number basis) the steady state 
response in orders of e: 



approx 

Po 



2c 



1 - 0(e 2 



\ 



-2A+17 1 

8e 3 (4 X +2A-z 7 ) 

(2A-i7) 2 (2x+2A-i7)(2A+i7) 

0(e 3 ) 



0(e 2 



(2A-i7)(2x+2A-i7)(2A+i7) 

0(e 4 ) • 
0(e 4 ) 0(e 5 ) 0(e 6 ) / 



(32) 
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FIG. 8: Response |(o)| for low damping, showing the emergence of the line at zero detuning. Here 7 = .01, \ = 1> an d from 
bottom to top, e = .002, .003, .004, .005, .006, .007, .008, .009, .01. At the low excitation level the response is Lorentzian and 
corresponds exactly to the linear response of a qubit; the higher lying levels play no role in the response. The hole in the center 
of the line at the highest e is also a qubit effect, showing the incipient saturation of the qubit response. The asymmetry of the 
line around zero detuning occurs because of participation of higher levels, and would be the first indication that the system is 
not a qubit as the excitation level is raised. 
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FIG. 9: Response |(a)| for low damping, showing the emergence of the line at detuning A = — x- Here 7 = .01, x = 1; and, 
from bottom to top, e = .004, .006, .008, .010, .012, .014, .016, .018, .02. The structure emerges as a Fano resonance out of the 
continuum of the far off resonant response at zero detuning. 



The matrix is Hcrmitian, and we have omitted the entries above the diagonal. We have only shown in detail the 
matrix elements (1|pq PPIOX |0) and (2|/5Q Pprox |l), since these are the only ones that contribute to the the observable (a); 
but, as we note, there are many non-zero unobservable components to this steady-state density operator. 

We also note the expression for the total response that can be extracted either from this or via a Taylor expansion 



ofEq. (20) 



2£ 32 * £3 + 0(e«). (33) 



w 2A-i7 (2A-i7) 2 (2x + 2A-i7)(2A + i7) 

First we turn our attention to the resonant structure near zero detuning A = 0, shown in Fig. [8| For small e, as 
the expression in Eq. (33) shows, the response is a simple Lorentzian, just as expected for a driven two-level system. 
In this regime, the oscillator behaves exactly like a qubit, which would, in the transient regime, respond as any two 
level system would to Rabi excitation or to Ramsey or spin-echo protocols. When e become of order 7, new structure 
appears. The top of the Lorentzian peak becomes flat and then exhibits a growing hole. One might say that at this 
point the higher levels of the oscillator are participating, but in fact this is very close to the behavior shown by a 
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two level system in the regime of spectroscopy probed at high power. As recently shown in [T7], saturation of the 
two level system causes just this same evolution into a double peaked structure. But some features of this regime are 
indicative of the participation of higher levels. The development of line asymmetry is definitely ascribed to higher 
levels; this asymmetry diminishes in the limit of large x- 

Mathematically, we see from Eq. (32) that these new structures in the zero-detuning line all arise from the 0(e 3 ) 



contributions. Note that only odd orders contribute observably to (a), consistent with the fact that the response 
should be an odd function of e. The 0(e 3 ) response involves all of the first four levels, so in some sense the line at this 
excitation level is indicative of four-level system response. But it is not clear that one need adopt this point of view, 
since the 0-3 matrix element of the steady state response is not visible to our problem, although one could construct 
a more complex observable that is sensitive to it. 

We next turn to an examination of the resonant structure near A = — %, at the "1 — > 2 transition" . The linear 
response has no resonant structure at this detuning, but the e 3 part of the response is resonant here. But the linear 
response has a uniform far off resonant background at this detuning, and resonance emerges our of this background 
(see Fig. [9|. Thus, the behavior here is very different in detail from the no-detuning line, although mathematically, 



they arise from two different poles of the same matrix elements in Eq. ( 32 ) . 

Furthermore, it is easy to show that because of interference of the resonant response with the background, the line 
that emerges in Fig. [9] has a Fano line shape [TS] 

(x - q) 2 

| (a) | ~ | (a) | background + C ^ — ^ (34) 

with suitably normalized frequency x. The Fano parameter q is easy to calculate: 

- 2 X 



V 2 X 2 + 7 2 - 7 



(35) 



Since we are in a regime where 7 << x, this Fano parameter is close to — 1. q = — 1 is among the most asymmetrical 
of the Fano line shapes [TU], in striking contrast to the initial symmetry of the distorted zero-detuning line. This Fano 
line exhibits a trough on the high-frequency side of the line. This is the beginning of the trough that forms a scallop 
from one "stalactite" to the next, and so it is the structure that has continuously evolved from the semiclassical trough 
that occurs at high damping. It is curious that the classical interference arising from intermittency at high damping 
evolves continuously into a Fano destructive interference, which we think of as entirely quantum mechanical. 

We will say much less about the resonant structures that arise at higher detuning, except to make a couple of 
observations: All higher lines emerge in a similar way to the A = — \ line, with a Fano line shape with the same 
asymmetry (so, evidently, a similar q, which we have not calculated). Each emerges first at successively higher order 
in the response; so, the resonant line at A = — n\ arises first at order 0(e 2n+1 ). The calculation shows that the first 
appearance of the line from the background (as measured by the first appearance of a zero slope near the resonance 
frequency) occurs at an excitation level that goes like 



^onset 



= f(n)x 1 -"l", (36) 



with some increasing function f(n). We note that the dominant background is always the off resonant response of 
the zero-detuning line, rather than the background of any nearer by lines, no matter how large the value of n. So, 



from Eq. (361 we see that it is only the n — and n = 1 lines that first show nonlinear structure at e ~ 7; structures 



at all higher detunings occur only at higher power {e.g., for n = 2 the onset is at y/xy). 



VI. DISCUSSION 



It is clear that the present study is only a piece of the full picture in the theoretical modeling of the complete 
input-output relations for the excitation of qubit structures. The potential complexity of the decohering environment 
of these structures has barely been touched on here, and remains a frontier area of research. For actual applications, 
more of the essential further work will be in the time domain rather than the frequency domain; pulsed excitation of 
qubits, both for gate operations and for measurements, will definitely be the study of more extensive modeling as this 
experimental art develops further. 

Still, the present simple study has revealed an intricate complex of phenomena. One example of a further study 
that it possible as a direct follow on to this work is the more detailed examination of the validity and limits of the 
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rotating wave approximation. The simplicity of the mathematics here should permit a definite if limited answer to 
questions about this. These questions often do not have such a clear answer in a more general setting. 

Another question raised by the present study concerns the emergence of classical bistability. As we have seen, the 
quantum theory says that there is a continuum of states (convex combinations of p+ and pJ) that are metastable. 
Which two of them does the classical physics pick out? We have offered a conjecture here that the ones of locally 
minimum von Neumann entropy (which happen to be p + and p_ themselves) are the classical states. Their low 
entropy is associated (albeit indirectly) with their high degree of localization of the Wigner function in the complex 
a plane. This in turn could make them candidates as "pointer states." That is, it is possible that the right view of 
the quantum environment would make these states emerge naturally. Thus, even the most fundamental questions can 
have a profitable launching point from our simple model calculations. 
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